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EFFICIENT  MODEL  FOR  HF  LASERS  WITH 


ROTATIONAL  NON  EQUILIBRIUM 


Rotational  relaxation  and  rotational  hole-burning  play  an  important 
role  in  determining  the  degree  of  saturation  and  the  rate  at  which  energy 
can  be  extracted  from  a chemical  - laser  medium.  Model  predictions  of  laser 
output  energy,  spectral  content,  and  performance  under  line  - selected 
operation  are  more  realistic  when  rotational- relaxation  mechanisms  are 
included.  Mode  structure  and  mode  competition  within  an  optical  cavity  are 
directly  related  to  the  distribution  of  gain  and  losses  within  the  cavity. 
Laser-output  beam-quality  predictions  should  be  made  with  a model  that  can 
accurately  predict  the  gain  distribution  within  the  active  medium.  Satisfacto- 
ry prediction  of  gain  distribution  requires  that  the  effects  of  rotational  non- 
equilibrium be  considered. 

Rotational  relaxation  has  been  modeled  either  by  means  of  detailed 

16  6 T 

rotational  kinetics  ’ or  by  use  of  a rotational  relaxation  time  constant.  ’ 

The  detailed  kinetics  model  is  more  realistic,  and  appears  to  have  the 

potential  of  being  able  to  predict  time  histories  of  the  laser  transitions  in 

all  the  comprehensive  detail  that  is  seen  in  experiments.  In  contrast,  a time 

constant  model  seems  to  display  the  dominant  features  of  this  time-resolved 

spectra  while  missing  some  fine  details.  Preliminary  comparisons  of  laser 

output  energy  and  spectral  energy  distributions,  however,  have  not  shown 

substantial  differences  in  the  two  modeling  approaches.^  More  exhaustive 

comparisons  between  the  models  and  with  experimental  spectra  need  to  be 


carried  out  in  order  to  identify  regimes  where  the  additional  realism 
afforded  by  the  detailed  kinetics  model  will  be  required. 

A difficulty  with  HF  laser  models  treating  finite  rate  rotational 
relaxation  is  their  need  for  large  amounts  of  computation  time.  The  method 
described  here  increases  the  efficiency  of  models  in  which  the  time-constant 
approach  is  used.  This  new  scheme  drastically  reduces  the  computational 
time,  with  little  loss  of  accuracy.  In  a comparison  with  a comprehensive 
earlier  model  (ROLAX)^  for  chain- reaction  HF  lasers,  the  number  of 
integration  variables  can  be  reduced  from  179  to  65,  and  computation  times 
are  typically  shorter  by  a factor  of  10  or  more.  For  the  two  models,  pre- 
dictions of  pulse  energy,  pulse  length,  and  spectral  energy  distribution 
typically  agree  within  5%. 

Formulation  of  the  new  model  through  our  Eq.  (4)  is  identical  to  that 
given  in  Ref.  7.  This  material  is  recapitulated  very  briefly.  The  major 
assumptions  are:  (1)  Kinetic  processes  are  represented  by  164  reactions, 

g 

based  on  the  survey  by  Cohen.  (2)  The  reacting  mixture  is  homogeneous 
and  contained  in  a Fabry- Perot  laser  cavity.  (3)  Lasing  is  assumed  to  be 
always  in  the  P-branch  and  restricted  to  the  first  fifteen  transitions  in 
each  of  the  lowest  six  vibrational  bands.  (4)  All  transitions  are  assumed 
to  have  low  initial  intensities,  which  grow  if  the  gains  rise  above  threshold. 

HF  vibration- rotation  state  populations  N(v,J)  are  computed  from  the 
kinetic  reactions,  with  the  effects  of  lasing,  absorption,  and  spontaneous 
emission  taken  into  account. 
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The  gain  coefficient  is 


hN 


a(v,  J)  = «c(v,J)  <pc  B(v,  J) 


gl 

— N(v  + 1,  J - 1)  - N(v,  J) 
gu 


(1) 


where  v and  J are  quantum  numbers  of  the  lower  level  of  the  transition, 

u)  (v,  J)  is  the  wavenumber  of  the  transition,  B(v,  J)  is  an  Einstein  coefficient, 
c 

and  g|  are  upper  and  lower  state  degeneracies,  respectively,  and  N(v,J) 
are  state  populations.  Both  Doppler  and  pressure  broadening  are  considered 
through  use  of  a Voigt  profile  <p which  is  evaluated  at  line  center. 

The  rate  of  photon  extraction  per  unit  volume  is 

*rad(v>J)  = a(v'J)  f(v'J)  (2) 

where  f(v,  J)  is  photon  flux.  The  rate  equation  for  f is 

= c(t)  [a(v'J)  - ath]  f(v*J)  {3) 

where  c is  the  speed  of  light  and 

“th  = " ~ZU  ln  (rorl)  (4) 

where  L is  the  length  of  the  active  medium,  £ is  the  mirror  spacing,  and 

R_  and  Rx  are  mirror  reflectivities. 

U .Li 
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For  brevity  and  clarity,  the  new  rotational  relaxation  model  is  * 

explained  in  terms  of  a single  vibrational  band  laser.  A discussion  of  the 
formulation  for  a multiband  laser  with  cascading  will  be  given  elsewhere 
because  of  its  length.  Numerical  results  given  here,  however,  are  for  the 
full  multiband  model. 

Let  N^d)  and  (J1)  denote  the  upper  and  lower  state  populations, 
respectively,  of  a laser  transition  with  photon  flux  f.  The  rate  equations  for 

the  populations  may  be  written  ‘ 

dN  (J)  Ne(J)  - N (J) 

— a* af  + — nn — (5a) 

u 

« 

and 

dN  (J')  Nf(J')  - N#.(J') 

it — = a{  + fyrn — (5b) 

where  N ^(J)  is  the  instantaneous  rotationally  equilibrated  concentration  of 

N.  (J)  and  t.(J)  is  the  time  constant  for  rotational  relaxation  of  level  N.(J). 

11  1 

The  rotational  relaxation  mechanism  represented  by  Eq.  (5)  is  ident- 
ical to  the  formulation  used  in  Ref.  7,  except  that  preferential  (i.  e.  , J- 
selective)  pumping  and  deactivation  are  ignored.  Here,  individual  rotational 
states  are  fed  (or  depleted)  only  through  rotational  relaxation  or  lasing. 

Pumping  and  vibrational- state  deactivation  are  manifested  indirectly  through 
their  effect  on  N°  terms. 
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The  relaxation  time  t is  assumed  to  depend  on  the  energy  separation 

7 

of  adjacent  rotational  states.  The  form 


Tj~  Aj  exp 


is  used,  where  Aj"  is  proportional  to  the  frequency  for  binary  collisions  with 
all  partners,  B is  a constant,  and  AEj  is  the  energy  separation  of  adjacent 
rotational  states. 

For  P-branch  transitions  a useful  approximation  is 


T = T . = T 

U l 


For  this  single-band  system,  Eq.  (1)  may  be  written 


S, 

a(J)  = <p  — — N (J)  - N|  (J‘) 

s u 


where  represents  all  quantities  outside  the  brackets  in  Eq.  (1).  A rota- 


0 

tionally  equilibrated  gain  a is  defined  as 


'g 

ae(J)  = <p  / N®(J)  - Nj(J') 
8u 


Differentiating  Eq.  (8)  with  respect  to  time,  and  substituting  Eqs.  (5) 


and  (7)  through  (9),  one  obtains 


= -<•(—  + i)  at  + 0,18 

dt  V®„  / 


’ ''  - ' 


At  saturation,  a ~ = constant.  Therefore,  da/dt  = 0 and  Eq.  (10)  can  be 

solved  for  a to  obtain 


where  the  "saturation  flux"  f is  defined  by 

s ’ 


(ID 


(12) 


Although  Eq.  (11)  was  derived  with  the  use  of  the  saturation  assumption,  it 

0 

is  clearly  valid  for  the  case  f « f since  it  reduces  to  a = a . The  calcula- 

e 

tions  in  Ref.  7 from  a rotational  equilibrium  model  (a  = a ) indicate  that  the 
saturation  process  is  rapid  enough  such  that,  prior  to  saturation,  the  condi- 
tion f«  f generally  holds.  Thus,  Eq.  (11)  may  be  taken  to  apply  at  all 
times  with  very  little  error.  This  result  makes  it  possible  to  predict  the 
nonequilibrium  gains  of  the  individual  vibrational-rotational  transitions 
directly  from  the  vibrational  populations  and  the  photon  fluxes  of  the  corre- 
sponding transitions.  The  need  to  compute  instantaneous  rotational  state 
populations  is  avoided,  thus  eliminating  them  as  required  integration 
variables. 

A rotational  bottle-necking  effect  is  clearly  demonstrated  by  Eq.  (11). 

The  observed  laser  intensity  depends  on  the  magnitude  of  the  saturation 

intensity  f , i.e.  , for  small  f , f must  remain  small.  Since  f is 
s s s 
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proportional  to  1/t,  the  rate  at  which  energy  is  pumped  into  (or  extracted 
from)  the  laser  inversion  is  limited  by  rotational  relaxation. 

From  Eq.  (6),  t increases  with  increasing  J.  Therefore,  f is  a 
J s 

decreasing  function  of  the  rotational  quantum  number.  Rotational  bottle- 
necking becomes  stronger  as  J increases.  At  the  same  time,  because  of 
the  increased  rotational  energy  spacings,  collisional  quenching  is  believed 
to  be  less  rapid  with  increasing  J.  An  energy  trapping- like  characteristic 
results,  which  may  be  a contributing  factor  in  the  comparatively  long  but 
low-powered  laser  emissions  observed  for  high  J transitions? 

The  preceding  formulation  was  extended  to  a multilevel  system,  and 
its  predictions  compared  with  those  of  the  comprehensive  rotational  non- 
equilibrium model  of  Ref.  6.  The  major  distinction  between  the  two  models 
is  in  the  calculation  of  gain.  The  model  described  in  Ref.  6 keeps  track  of  the 
population  of  each  v,  J state  and  uses  these  in  Eq.  (1)  to  determine  the  gain 
for  each  transition.  The  present  model  keeps  track  of  only  the  vibrational 

state  populations.  The  Boltzmann  distribution  then  is  applied  to  obtain  the 

0 

equilibrated  gains  a , which  are  subsequently  utilized  in  the  multiband 
counterpart  of  Eq.  (11)  to  obtain  the  nonequilibrium  transition  gains.  The 
kinetics  utilized  in  both  models  are  identical.  For  consistency,  the  relaxa- 
tion time-constant  version  of  the  Ref.  6 model  is  used. 

A comparison  of  the  predicted  pulse  outputs  of  the  two  models  for  an 
atmospheric  pressure  H^-F^  laser  is  given  in  Fig.  1.  The  pulses  are  nearly 
identical,  even  to  the  details  of  the  individual  "bumps"  in  the  pulse  shapes. 


J.  S.  Whittier  and  R.  Hofland,  private  communication,  The  Aerospace 
Corporation,  El  Segundo,  California. 


In  Fig.  2,  the  model  predictions  of  time  resolved  spectra  for  the  same  laser 
case  are  compared.  The  horizontal  bars  denote  the  intervals  during  which 
lasing  is  observed  for  the  transitions  designated  on  the  vertical  scale,  where 
J is  the  rotational  quantum  number  of  the  lower  level  of  the  transition. 

Again,  the  predictions  are  nearly  identical.  This  consistency  is  also  observed 
in  the  relative  energy  distributions  of  the  individual  vibrational- rotational 
transitions  of  the  six  lasing  bands.  The  amount  of  information  attainable 
by  the  use  of  either  model  is  the  same.  As  demonstrated  in  the  figures,  the 
predictions  are  very  close.  The  model  of  Ref.  6 explicitly  computes  the 
populations  of  the  individual  rotational  states.  In  the  present  model,  this 
information  can  be  inferred  from  the  calculated  transition  gains  and  vibra- 
tional populations. 

The  cost  savings  achieved  by  the  present  model  are  considerable.  For 

the  case  described  here,  the  execution  time  for  the  Ref.  6 model  was  1640 

sec  on  the  CDC  7600  computer.  For  the  present  model,  the  time  required 

is  just  132  sec.  The  execution  time  of  the  corresponding  model  without 

7 

rotational  relaxation  features  is  104  sec.  Additional  reductions  in  the  pre- 
sent model  run  time  are  anticipated  as  the  study  progresses. 

Since  the  rotational  nonequilibrium  model  utilizing  the  technique 
described  herein  costs  only  slightly  more  to  run  than  the  corresponding 
equilibrium  model,  laser  calculations  with  rotational  relaxation  taken  into 
account  can  now  be  carried  out  routinely.  Furthermore,  parametric 
studies  can  be  performed  that  will  improve  understanding  of  the  importance 
of  rotational  relaxation  in  laser  predictions.  The  technique  is  applicable  to 
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Comparison  of  Time-Resolved  Spectra 
of  Present  Model  and  ROLAX 


both  pulsed  and  cw  laser  calculations.  Of  particular  interest  is  its  possible 
application  to  wave  optics  modeling  of  cw  laser  beam  quality,  where  the 
added  spectral  realism  afforded  by  considering  rotational  nonequilibrium  can 
now  be  obtained  at  a much  more  reasonable  cost. 
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developing  space  and  missile  systems.  Expertise  in  the  latest  scientific  devel- 
opments is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aerophysics  Laboratory:  Launch  and  reentry  aerodynamics,  heat  trans- 
fer, reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
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pheric  optics,  cnemical  reactions  in  polluted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo- 
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Electronics  Research  Laboratory:  Electromagnetic  theory,  devices,  and 
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and  airglow;  magnetospheric  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  magnetosphere;  solar  physics,  studies  of  solar  magnetic 
fields;  space  astronomy,  x-ray  astronomy;  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
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